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The NA60 experiment has studied low-mass muon pair production in In-In collisions at 158 
AGeV with unprecedented precision. Thus, there is hope that the in-medium modifications of the 
vector meson spectral function can be constrained more thoroughly than from previous experiments. 
Towards this goal, we present a study comparing different models for the in-medium modification 
of vector mesons. In particular, we investigate what can be learned about collisional broadening 
by a medium at finite temperature and baryon density and what model independent constraints 
can be obtained from the data. Of special relevance is a comparison to py-cuts of the dilepton 
spectra. Those not only provide independent tests of the evolution model and p-r-dependence of the 
spectral function but also unambiguously decide if the observed yield in the higher mass spectrum 
M > 1 GeV is mainly due to the hot phase of the evolution or four pion annihilation processes. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

The properties of vector mesons, especially the p-meson, 
are predicted to change in the presence of a hot and dense 
nuclear medium. The CERES and HELIOS-3 collabo- 
rations demonstrated that central nucleus- nucleus colli- 
sions result in a strong enhancement of the emission of 
low-mass dileptons as compared to scaled proton-nucleus 
and proton-proton collisions 0, S B ■ 
One explanation of these experimental findings is based 
on broadening of the p-meson due to interactions with 
the surrounding medium. Considerable broadening due 
to the mesonic and the baryonic component as a function 
of temperature and density has been found in various 
hadronic model calculation, for reviews see e.g 0, 0, 0- 
We will focus here on three specific calculations, namely 
the application of thermal field theory in a chiral frame- 
work a mean field model for hot and dense nuclear 
matter |9j and a non-perturbativc ^-functional approach 

The mass resolution of the dilepton spectra of the CERES 
experiment was not high enough to conclusively decide 
which of these scenarios is realized in nature. Unfortu- 
nately, not even the seemingly straightforward question 
whether a mass shift of the p-meson due to the Brown- 
Rho scaling law is realized in nature or not could be an- 
swered conclusively 0. 

The high precision di-muon data of the NA60 experiment 
now provide a possibility to put strong constraints on 
models of the in-medium effects on vector mesons. To- 
wards this aim, we conduct a comparison of the different 
model calculations with the data. 
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Our strategy is as follows: Based on an evolution model 
successfully used to describe a plethora of observables in 
Pb-Pb collisions [l^, we rely on scaling arguments to 
make a sophisticated attempt at describing the evolu- 
tion of matter in In-In collisions. Ideally, one would need 
to tune this evolution scenario to hadronic observables, 
however in practice these are not yet sufficiently well de- 
termined experimentally. We will address an estimate 
of the error introduced by our lack of knowledge of the 
underlying evolution later on. 

Using this evolution model, we then calculate the re- 
sulting dilepton spectrum for different assumptions of 
in-medium vector meson properties. A first test of the 
evolution model is the comparison with the integrated pr 
data. We also show results for low and high px cuts. To 
the extent that they are reflected in the data, they will 
lend a posteriori further confidence to the scaling laws 
used to determine the evolution. 

We argue that different regions of pt in the di-muon spec- 
tra probe different physics, because the px dependence of 
the observed spectrum is in a characteristic way sensitive 
to the evolution of the flow profile and temperature. 
Studying different regions of px also provides the pos- 
sibility to decide if the observed excess in the mass re- 
gion (0.9 < M < 1.5 GeV) is mainl y d ue the emission 
from a hadronic gas (as suggested in in later stages 
at lower temperatures (accompanied by higher flow) or 
from early stage radiation from a hot medium (accom- 
panied by lower flow). Our calculation of the medium 
evolution suggests that those contribution are predomi- 
nantly due to the early times of the medium evolution 
whereas the authors in [T^ | argue that four pion anni- 
hilation processes are required to account for the data. 
Since emission from this temperature region is greatly ef- 
fected by flow, considerable differences in the predictions 
of the pT-dependence of the spectra in both scenarios are 
expected. 
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The paper is organized as follows: In the next section 
we discuss how the dilepton spectrum is calculated. We 
especially explain how we employ scaling arguments to 
readjust a model that was able to describe hadronic 
transverse mass spectra and HBT correlation measure- 
ments and dileptons in 158 AGeV Pb-Pb and Pb-Au col- 
lisions at SPS to the In-In system. We then briefly dis- 
cuss the calculations that were employed to descripe the 
spectral function in the regime T < T c and T > T c , 
respectively. We proceed to explain the origin of p T de- 
pendence and the technique employed to account for flow 
effects on the vacuum p signal. We present the results of 
our model calculation and furthermore discuss test cal- 
culations that illustrate the role of the evolution model. 
Ther we also present the physical mechanism mainly re- 
sponsible for the high mass (M ss 0.9 — 1.5 MeV) excess. 
We conclude the paper with a general discussion of our 
results. 



II. CALCULATING THE DIMUON SPECTRUM 

The emitted spectrum of dileptons from an evolving ther- 
malized system 01 can be found from the convolution, 
schematically: 



are implemented in a more than just schematic way in 
order to allow for a meaningful comparison of theoretical 
spectral functions with the data. 

Every measurement of the dilepton spectrum always 
probes a combination of evolution and spectral function 
as apparent from Eq. It is therefore necessary to de- 
termine the evolution as reliably as possible in order to 
conclusively address the question of the behavior of the 
in- medium spectral function R(q, T, ps)- We now discuss 
each of the ingredients in turn. 



A. The evolution model 

A detailed description of the model used to describe the 
fireball evolution is found in [l2T |. The main assumption 
for the model is that an equilibrated system is formed a 
short time ro after the onset of the collision. Further- 
more, we assume that this thermal fireball subsequently 
expands iscntropically until the mean free path of par- 
ticles exceeds (at a timescale r/) the dimensions of the 
system and particles move without significant interaction 
to the detector. 

For the entropy density at a given proper time we make 
the ansatz 



d 3 N 



dMdpTdrj 



evolution ® 



dN 
d 4 xd 4 q 



acceptance , (1) 



where the rate calculated from the averaged virtual pho- 
ton spectral function R(q,T, ps) as 



s(r,r) s ,r) = NR(r,r) ■ H( Vs ,t) 



(3) 



with t the proper time as measured in a frame co-moving 
with a given volume element and i?(r, t), H(t] s ,t) two 
functions describing the shape of the distribution and N a 
normalization factor. We use Woods-Saxon distributions 



dN a 2 R(g,T, PB ) 
d 4 xd 4 q 12tt 4 eAv«" - 1 ' U 

Here, q represent the four-momentum of the emitted 
muon pair, T is the temperature of the emitting volume 
element and pb its baryon density. The fireball evolution 
encodes information on the evolution of radiating vol- 
ume, temperature T and baryon chemical potential pb, 
local transverse flow velocity vt, longitudinal rapidity r\ 
and in the late evolution stage chemical non-equilibrium 
properties such as pion chemical potential p^ and Kaon 
chemical potential p,K- 

The acceptance is highly non-trivial in the the case of 
NA60. Therefore for all instances where we compare our 
calculations with the data, a full Monte Carlo simula- 
tion of the detector has been carried out by the NA60 
collaboration with our calculations as input. 
Although a Lorentz-invariant quantity like the invariant 
mass spectrum dN/dM is in principle independent from 
the flow profile of the medium, the resulting invariant 
mass spectrum after acceptance folding is not: The ac- 
ceptance has pT-dependence and thus the hardening of 
the momentum spectra due to flow alters the shape of 
the invariant mass distribution. Therefore it is manda- 
tory that both the medium evolution and the acceptance 



R{r,r) = lJ fl + cxp V / c(r) 1) 

H(rj s ,r) = l/ (l + e,J^m). 

to describe the shapes for a given r. Thus, the ingredi- 
ents of the model are the skin thickness parameters d ws 
and rj ws and the parametrizations of the expansion of the 
spatial extensions R c (t), H c (t) as a function of proper 
time. For a radially non-relativistic expansion and con- 
stant acceleration we find R c {t) = Ro + ^t 2 . H c (t) is 
obtained by integrating forward in r a trajectory origi- 
nating from the collision center which is characterized by 
a rapidity t] c (t) = T]o + a v r with i] c = atanh v c z where 
v c z is the longitudinal velocity of that trajectory. Since 
the relation between proper time as measured in the co- 
moving frame and lab time is determined by the rapidity 
at a given time, the resulting inte gral is in general non- 
trivial and solved numerically (see [l3| for details) . Ro is 
determined in overlap calculations using Glauber theory. 
rjo is the initial size of the rapidity interval occupied by 
the fireball matter, and a± are free parameters, but 
we choose to use the transverse velocity = a±Tf and 
rapidity at decoupling proper time rj' — rjo + a^Tt as 
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parameters instead. Thus, specifying r)o,r)f,Vj_ and r/ 
sets the essential scales of the spacetime evolution and 
d ws and rj ws specify the detailed distribution of entropy 
density. 

For transverse flow we assume a linear relation between 
radius r and transverse rapidity p = atanh vj_(t) = 
r/R c (r) ■ p c (r) with p c (r) = atanh aj_r. 
We allow for the possibility of accelerated longitudinal 
expansion which in general implies r\ ^ r\ s |l2j. Here, 
r? = i In Po+p - denotes the longitudinal momentum ra- 
pidity of a given volume element. We can parametrize 
this mismatch between spacetime and momentum rapid- 
ity as a local A?; = ?/ — t? s which is a function of r and 

Vs- 

The model parameters have been adjusted to hadronic 
transverse mass spectra and HBT correlation measure- 
ments in 158 AGeV Pb-Pb and Pb-Au collisions at SPS. 
The framework was then used to successfully describe 
photon and dilcpton emission and charmonium suppres- 
sion Since at this point no freeze-out analysis or 
HBT correlation data are available from NA60 for the 
In-In system, our strategy is to start with the parame- 
ter set derived for Pb-Pb collisions and use geometrical 
scaling arguments to go to In-In. 

We infer the change in total entropy production from 
the ratio of charged particle rapidity densities dN c h/di] 
at 30% peripheral Pb-Au collisions with 2.1 < rj < 
2.55 measured by CERES 0, H, and semi-central 
In-In collisions at 77 = 3.8 measured by NA60 [Tlj . 

dN ch~ In / dN ch~ Pb = °- 68 - Thc number of participant 
baryons and the effective initial radius (for the sake of 
simplicity, we map the overlap area in non-central colli- 
sions to a circle with the same area) are obtained with nu- 
clear overlap calculations. We assume that the stopping 
power (which determines the width of the inital distribu- 
tion of entropy and baryon number in rapidity) scales ap- 
proximately with the number of binary collision per par- 
ticipant, Nbin/Npart- This determines the scale t]q. The 
shape parameters d ws and rj ws have no great impact on 
electromagnetic emission into the midrapidity slice, they 
primarily govern the ratio of surface to volume emission 
of hadrons. Thus, we leave them unchanged from their 
value determined in [l2T | . Under the assumption that the 
physics leading to equilibration is primarily a function of 
incident energy we keep the formation time r = 1 fm/c 
as in Pb-Pb collisions. We stress that the final results 
exhibit no great sensitivity to the choice of either 770 or 
To except in the high M, high px limit. 
The biggest uncertaintly is the choice of the decoupling 
temperature Tp, Due to the fact that the In-In system 
is smaller than Pb-Pb, we expect a higher decoupling 
temperature. However, without simultaneous knowledge 
of HBT correlations and transverse mass spectra an un- 
ambiguous answer cannot be obtained. We tentatively 
choose Tp = 130 MeV for the time being. Our lack of 
detailed knowledge translates primarily into an uncer- 
tainty in the normalization of thc final results (which de- 
pend on the length of emission duration which in turn is 



increased for a lower decoupling temperature). However 
we find that our model describes the absolute normal- 
ization rather well with this T/, so we do not perform 
further adjustments. 

We fix the flow velocity at Tc by choosing the value of 
transverse flow obtained at T = 130 McV in our Pb-Pb 
model in the same centrality class. There is likewise an 
uncertainty (which is to some degree correlated with the 
decoupling temperature) since flow determines the slope 
of hadronic (and to a lesser degree electromagnetic) pr 
spectra. 

The equation of state in the hadronic phase and the off- 
equilibrium parameters p^PK are inferred from statis- 
tical model calculations as described in 0, 0| . The re- 
sulting fireball evolution is characterized by a peak tem- 
perature of about 250 MeV, a lifetime of ~ 7.5 fm/c and 
a top transverse flow velocity of 0.55c at decoupling. The 
scales are quite similar to those obtained in a hydrody- 
namical evolution model in |l6l | but differ considerably 
from [13 . 



B. Hadronic model calculations 

The effective degrees of freedom below the phase transi- 
tion temperature are color singlet bound states, namely 
hadrons. Thc photon couples to the lowest-lying dipole 
excitations of the vacuum, thc hadronic J p =l~ states: 
the p, uj and <f> mesons and multi-pionic states with thc 
appropriate quantum number. The dilepton data from 
CERN-SPS has been the driving force for theoretical at- 
tempts to understand how vector-meson properties are 
modified in hot and dense strongly interacting matter. 
Especially the p-meson has been the subject of many 
theoretical investigations due to its prominent role for 
dilepton emission. 

Different effective hadronic models and techniques for 
calculations of the properties of hadronic matter near 
and below the phase boundaries have been suggested. 
For reviews of thc intense theoretical activities, see e.g. 

ilQ 

Most of those calculations constrained by phenomenolog- 
ical data such as the hadronic and electromagnetic decay 
widths of the particles predict substantial broadening in 
matter with comparably small shifts of thc in-medium 
masses. 

We study in the present publication the implications of 
three different models for the calculation of the observed 
dilcpton spectra and compare their predictions for the 
changes of thc in-mcdium properties of thc vector mesons 
quantitatively with the NA 60 data below. 

The first model is based on 0, where the electro- 
magnetic current-current correlator has been computed 
in the so called improved vector meson dominance model 
combined with the chiral dynamics of pions and kaons. 
The irreducible photon self-energy in this model can be 
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directly related to the vector meson self-energy: 



imn 



F V {q) 
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(5a) 
(5b) 



here my are the in-medium renormalizcd vector meson 
masses, g v is the photon vector meson coupling, and g 
the coupling of the vector meson to the pseudoscalar 
Goldstone bosons. This equation is valid in the model 
for vanishing three momentum q and the calculation of 
the electromagnetic current-current correlator is applied 
in this limit. For the evaluation of finite baryon den- 
sity effects which are relevant at SPS conditions the re- 
sults from Klingl, Kaiser, and Weise ^| are used. Ther- 
mal broadening of the p, 0, and uj were calculated by 
Schneider and Weise as modifications on top of the finite 
baryon density effects using perturbativc methods |sl ll7j . 
In these calculations it is assumed that the temperature 
and density dependences of the vector meson self-energy 
factorize. This amounts to neglecting contributions from 
pion-nucleon scattering where the pion comes from the 
heat baths. Formally speaking matrix element contribu- 
tions such as (TrN\Tj^(x)j fJ, (x)\TrN) are not taken into 
account. 

A second approach Q is based on a different idea: The 
method of thcrmoficld dynamics is used to calculate the 
state with minimum thermodynamic potential at finite 
temperature and density within the mean field sigma 
model with a quartic scalar self interaction. The temper- 
ature and density dependent baryon and sigma masses 
are calculated self-consistently. The medium modifica- 
tion to the masses of the u>- and p-mesons in hot nuclear 
matter including the quantum correction effects are then 
calculated in the rclativistic random phase approxima- 



tion. The decay widths for the mesons are calculated 
from the imaginary part of the self energy using the 
Cutkosky rule. The model also includes an additional 
p broadening contribution by scattering off by baryons. 

We also discuss a third model calculation where we 
solved truncated Schwinger-Dyson equations derived in 
a ^-functional approach incorporating a self-consistent 
resummation of the ir — p interaction in order to de- 
termine the thermal broadening of the p-meson, see 
|lOj and references therein. This ^-functional approach 
sclf-consistcntly takes into account the finite in-medium 
damping width of the pion in the thermal heat bath 
which contributes to the broadening of the p-meson 
[Tol fl9| . Effects of baryons are not yet included in a 
self-consistent approximation scheme. 

The solution in the ^-functional approach includes the 
full momentum dependence of the three-dimensional lon- 
gitudinal and transverse components of the p-mcson's 
spectral function which is especially important for a com- 
parison with the NA60 data measured at integrated and 
high p t . 

C. The partonic spectral function 



We describe the thermalized partonic evolution phase by 
means of a quasiparticle model. Details can be found in 
H, here we only summarize the essential findings. 

As long as the evolution is in a regime where the ther- 
modynamically active degrees of freedom are quarks and 
gluons, the timelike photon couples to the continuum 
of thermally excited qq states and subsequently converts 
into a charged lepton pair. 

The photon spectral function has been calculated at the 
one-loop level using standard thermal field theory meth- 
ods. The well-known leading-order result for bare quarks 
and gluons as degrees of freedom is: 




where q — (q° , q) is the four-momentum of the virtual 
photon, e/ the quark electric charge, fn the Fermi-Dirac 
distribution and my the quark mass of flavour /. This re- 
sult, however, holds only up to perturbative higher order 
corrections in g s that take into account collective plasma 
effects. 



Let us now assume that a quark quasiparticle couples 
to a photon in the same way as a bare quark (a form 
factor representing the 'cloud' of the quasiparticle could 
in principle also be included, but in absence of informa- 
tion about the detailed quasiparticle structure we ignore 
this point). For a gas of non-interacting quasiparticlcs, 
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the one-loop result for R is adequate, with input prop- 
erly adjusted. In the formalism outlined in all higher 
order QCD effects manifest in the thermal quasiparticlc 
masses ro(T), the thermal vacuum energy B(T) and the 
confinement factor C(T). Incorporation of the first two 
features in the calculation is straightforward. The bare 
quark masses in eq. © simply have to be replaced by the 
T-dcpendent quasiparticlc masses for each flavour. The 
thermal vacuum energy B(T) does not contribute to the 
dilcpton rate. 

The mechanism for dilepton production at tree-level is 
the annihilation of a qq pair into a virtual photon where 
the quark lines arc multiplied by the distributions /d(T), 
giving the probability of finding a quark or an antiquark 
in the hot medium. The incorporation of the confine- 
ment factor is hence straightforward: since it reduces the 
number of thermally active degrees of freedom, it also 
reduces the dilepton rate by a factor of C(T) 2 . 
Let us note at this point that QGP radiation is expected 
to be important only in the limit of dileptons at high 
M and high px- However, in this limit (which corre- 
sponds to the high temperature phase of the evolution) 
the quasiparticle description is not very different from the 
free quark description, thus we expect the overall uncer- 
tainty resulting from the quasiparticle treatment of the 
QGP to be small. 



D. The origin of px dependence 

The transverse momentum dependence of a given dimuon 
detected at midrapidity arises in this framework for two 
different reasons. First, the virtual photon created in 
a qq collision or the decay of a vector meson may not 
be purely timelike but have a finite 3-momentum (with 
respect to the thermal heat bath) instead. This induces 
a 3-momentum dependence into the in-medium spectral 
function, hence we have R(qo, q, T, pg) in general. 
However, flow effects alter the momentum dependence. 
While neither photons nor muons take part in flow them- 
selves, the matter they are emitted from is hadronic and 
participates in collective motion, hence the local rest- 
frame in which virtual photons are emitted is in gen- 
eral not the lab frame. Thus, (neglecting rapidity depen- 
dence for a moment) muon pairs emitted from a moving 
volume clement roughly gain an additional momentum 
Pflow ~ AfvT where M is the invariant mass of the pair 
and vt the local transverse flow velocity. 
For this reason, the experimentally measured pt does 
in principle not allow to make conclusions about the 3- 
momcntum of the virtual photon in the local rcstframc. 
This is an important point since in several model calcu- 
lations for the spectral function the simplifying assump- 
tion q = has been made. We have checked that this 
assumption leads to harder px spectra and overestimates 
the yield at high px significantly. In this case, only a 
comparison with the data in the range px < 0.5 GeV is 
meaningful. While this doesn't strictly guarantee px ~ 



in the local restframe, we may note that for small in- 
variant masses M < 0.5 GeV Pfl ow < 250 MeV which 
ensures a small momentum q for most of the integration 
over the spacetime integration region. We may also note 
that in the model masses above 1 GeV are dominated by 
radiation from hot regions, i.e. early times for which the 
transverse flow (and hence Pfl ow ) is small. 

We have checked within a scenario with full 3-momcntum 
dependence in the spectral function that the distortion 
of the lineshape due to neglecting q is acceptable (on the 
order of at most 20%) as long as the comparison with the 
data in the region px < 0.5 GeV is made, a comparison 
of these calculations with the full px integrated data is 
however not meaningful at all. 



E. The vacuum p-meson 

Additionally to the dilcpton radiation from the medium 
(which is dominated by the channel it + it~ —> p° — > 7* — » 
and its in-medium corrections) there is also radia- 
tion from electromagnetic decays from p mesons outside 
the medium after the thermalized evolution has ceased, 
i.e. after the hadronic freeze-out. 

In principle this radiation originates also from uj and <p 
meson vacuum decays, however, due to the small width 
of these states this contribution is experimentally very 
well identified and has been subtracted ■ The vacuum 
decay of p mesons is commonly referred to as 'cocktail p' 
and obtained by calculating the produced number of p 
mesons by a statistical model calculation. 

In such calculations, the number of p mesons is de- 
termined under the assumption that the relative abun- 
dances of hadron species are given by the thermal expec- 
tation value above a temperature Th (close to the tran- 
sition temperature Tq) but that subsequently only reso- 
nance decay (and regeneration) occur. In such a frame- 
work, the decay of heavy resonances leads to an over- 
population of pion phase space as compared to the pion 
phase space density in chemical equilibrium. This can 
be parametrized by the introduction of a pion chemical 
potential (i^ (which grows with decreasing temperature) 
and leads via the tightly coupled p irir channel to an 
effective p, p = 2p, n . 

However, while tracking of the resonance decays fixes the 
number of p-mesons at thermal decoupling, their momen- 
tum spectrum remains unknown. In particular, p mesons 
after chemical decoupling undergo not only resonance de- 
cay and regeneration but also participate in the collective 
expansion of matter. While this effect does not influence 
the integrated invariant mass spectrum, the acceptance 
leads to a significant influence of flow on the detected 
signal. 

Thus, we go beyond the statistical model formulation and 
compute the full momentum spectrum of p mesons after 
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decoupling using the Cooper-Frye formula 

E ^Pp~ = J2^J da » P ° XP [ T f J (7) 

with p 11 the momentum of the emitted p and defy an 
element of the freeze-out hypcrsurface (determined by 
the condition T = Tp) just as we would calculate the 
emission of any other hadron from the medium. We let 
this ensemble undergo electromagnetic decay with vac- 
uum width and lifetime. In this way, the vacuum p sig- 
nal is influenced by the fireball flow field at the time of 
decoupling as one would expect. 



III. RESULTS 

A. In-medium modifications from many-body 
calculations 

We present a comparison of the calculation with the data 
for pt integrated data using the spectral function ob- 
tained in the ^-functional approach for the 7r — p sub- 
system [Tol | in Fig. n Furthermore we present calcula- 
tions of different muon-pair p^-cuts, i. e. a high px- 
cut with pt > 1 GeV in Fig. ^ and a low p^-cut with 
p T < 0.5 GeV in Fig.0 

All calculations contain the vacuum p calculated using 
the Cooper-Frye formula as described in section III El 
Since most of the vacuum p contribution is generated 
during the final breakup of the fireball ('volume freeze 
out') the underlying flow is strong and shuffles strength 
to higher pt- This is the reason that the vacuum p con- 
tribution is almost negligible for pt < 0.5 GeV (compar- 
ison not shown), but it is important for integrated pt 
and pt > 1 GeV, see Fig. right panel. 
We note that the relative normalization of vacuum p and 
in- medium p appears to be approximately right. While 
the strength of the vacuum p parametrically scales like 
the freeze-out volume times the p lifetime V(Tf)xjj ac , the 
in-medium radiation scales like J V(r)d,T and is thus 
sensitive to the lifetime of the fireball. Thus the good 
description of the the intermediate mass region M w 
0.6 — 0.9 GeV is due to contributions of the vacuum and 
in-medium p and therefore a consistency test of the fire- 
ball model. This will be even more stringent for the high 
Pt cut data, because the relative contribution of the vac- 
uum p to the in-medium p is even more pronounced. 
We note that the agreement of the ^-derivable approach 
calculation with the data including the vacuum p is in 
very good agreement over the whole range in M. This is 
to some degree surprising and unexpected, as the model 
contains substantial amount of the non-perturbative dy- 
namics of the 7r — p interaction but is completely insensi- 
tive to the presence of other hadrons, especially baryons 
(which, based on perturbation theory, are expected to 
play a major role in the broadening of the p). 



One may take this as an indication that the generic shape 
of the spectral function and its 3-momentum dependence 
is right, not a proof that this is the main mechanism 
capable of producing such a spectral function or that the 
dynamics of baryons is unimportant. However, this result 
clearly indicates that it is necessary to investigate non- 
perturbative mechanisms for broadening in more detail. 
In Fig. [21 we also show the results expected from a chiral 
model calculation || and a thermal nuclear matter mean 
field theory calculation under the assumption q = 
which restricts the predictive power to the low pt region. 
The chiral calculation also contains an in-medium and 
u> contribution (as these have different mass/width in 
medium than in vacuum, such contributions would not be 
removed by the experimental subtraction of the vacuum 
4>, to.) One interesting finding here is with regard to the 
</i>-meson. Although in the pt integrated data there ap- 
pears some room for an contribution in-medium around 
the <f> mass, no such signal is apparent in the preliminary 
low pt data as presented by S. Damjanovic at Quark 
Matter 2005 [13]. The KEK-PS E325 collaboration has 
announced an observation of a considerable in-medium 
modification of the 0-mcson in 12 GeV proton-nucleus 
collisons that is compatible with the chiral model pre- 
dictions at low temperature and normal nuclear matter 
density 0, |2l| . Still one has to be careful: in-medium 
radiation through the p channel is primarily sensitive to 
the coupling of 7r + + ir~ <-> p whereas radiation from the 
4> probes K + +K~ <-> 4> and it is by no means self-evident 
that both need to be in thermal equilibrium at all times 
during the evolution. 



B. The role of the evolution model 

Looking at the pt integrated invariant mass spectrum, 
one is tempted to describe it as a vacuum p peak on 
top of an almost symmetrically broadened structure, and 
naively one would infer that therefore the in-medium 
p spectral function must show symmetrical broadening. 
However, it is rather easy to realize that this view has a 
major flaw: As Eq. J2J indicates, the emission rate is not 
given by the spectral function but by the spectral func- 
tion multiplied by a thermal distribution. Thus, if the 
spectral function exhibits symmetric broadening, the cor- 
responding rate cannot be symmetric for typical temper- 
atures reached in heavy ion collisions but falls off sharply 
with increasing M. Furthermore, in the actual calcula- 
tions, contributions from many volume elements with dif- 
ferent flow and temperature are sampled and, indeed, the 
calculated distribution d 2 N/dMdpT is anything but flat 
even for a symmetrically broadened spectral function. 
It is a lucky but merely accidental coincidence that for a 
certain class of evolution scenarios a flat spectral function 
after folding with the fireball evolution and the accep- 
tance returns a flat mass distribution over some range in 
M. Thus the data in this sense approximate the (T, pb)) 
averaged lincshapc of the spectral function. We illustrate 
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FIG. 1: Left panel: comparison of the NA60 dimuon data [Tl| with calculations within a self-consistent ^-functional approach 
for the 7r — p interaction [Tflj . Right panel: Calculation within a self-consistent ^-functional approach for the ty — p interaction 
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Q for q = 0, and a thermal nuclear matter calculation ('MR') @ for q = 0. Right panel: Comparing the evolution without 
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this property in Fig. [3 left panel with the scenario called 
'full evolution'. Here, we use the qq continuum spectral 
function (the strongest possible symmetrical broadening) 
throughout the whole evolution, i.e. above Tq and be- 
low Tc- As aparent from the figure, an approximately 
flat signal is recovered after the acceptance folding in the 
mass region 0.25 GeV < M < 1.2 GeV. 

However, this is not a general property of the accep- 
tance but demands a specific type of evolution. To il- 
lustrate this point, we integrate radiation from the evo- 
lution model only from regions where T < Tc (this could 
be pictured as the result of a long equilibration time To). 
In spite of the fact that the spectral function has the same 
continuum shape, the output is not flat in this case. It 
can clearly be observed that radiation from hot volume 



elements with T > Tc is increasingly important for the 
contribution to invariant masses above the p mass. 

This argument can be turned around: Even the maximal 
amount of broadening possible in a hadronic phase (i.e. 
the continuum spectral function) will not reproduce the 
data - if the normalization is such that the high M region 
is described correctly, the yield overshoots the invariant 
mass region below the p peak, if the low M region is 
described correctly, there is not enough yield above the 
p peak. We call this effect here therefore 'excess in the 
higher mass region' (M « 0.9 — 1.5 GeV). 

We illustrate this point further in Fig. [3J right panel 
where the importance of radiation from regions with 
T > Tc for the description of the invariant mass spec- 
trum above the p peak is clearly visible when the full 
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FIG. 3: Comparison of the NA60 dimuon data [Tl| for integrated pr with different assumptions about spectral function and 
model evolution. Left panel: assuming qq spectral function at integrated pr for the full evolution ('hot') and the evolution 
below Tc ('cold'). Right panel: Comparing the evolution below T c — 170MeV with the evolution for all temperatures. 



hadronic spectral function is used. 

This indicates that the most likely conclusion is that the 
data require radiation from a significant amount of mat- 
ter with T > Tc in order to account for the excess in 
the higher mass region. Our calculation of the medium 
evolution suggest therefore that those contributions are 
predominantly due to the early times of the medium evo- 
lution. 

The authors of argue differently and claim that a fur- 
ther medium effect, namely contributions from four pion 
annihilation processes, accounts for the observed high 
mass excess. Since contributions to the spectrum from 
four pion annihilation would be dominated by the part of 
the evolution where the pion-chemical potential is large 
(and the rate is augmented by the pion fugacity factor 
exp (4/ijr/T)), they have to assume that the process is im- 
portant at lower T w 120MeV. They furthermore need to 
have considerable chiral mixing of the free V and A cor- 
relators even at comparatively low temperatures. Their 
calculation estimates an upper limit on the in-medium 
four pion annihilation contribution to dilcpton emission 
|l3l 12 2i | . This upper limit of the four pion annihilation 
with chiral mixing contribution to the spectrum is barely 
enough to account for the high mass M > IGeV region 
of the pt integrated data in a medium evolution model 
with an initial temperature of 197 MeV. 

One might doubt on theoretical grounds that all of the 
assumption in favor of this upper limit of the four pion 
annihilation contribution are realized. But the question 
if the high mass excess is mainly due to emission at high 
temperatures (accompanied by insignificant flow) or low 
temperatures (accompanied by large flow) can be decided 
unambiguously experimentally. 

The important observation is that the px- dependence 
of the spectra at given M with flow or temperature is 
different. Flow simply pushes spectra to higher px and 
reduces strength at lower pt accordingly whereas tem- 



perature dominated spectra are exponential without re- 
ducing strength at low px- Pt spectra at large M profit 
more from high flow than low mass spectra, for temper- 
ature dependence the opposite is true. 
Given only a comparison with the px integrated data one 
can have similar spectra from both mechanism since one 
can trade high px contribution against low px contribu- 
tions by suitable model adjustments. The two different 
mechanism can therefore lead to comparable predictions 
in the region M > 1 GeV for all px spectra. 
This is different if one wants to describe low px and high 
px cuts at all M consistently with one mechanism. For 
this not only an inclusion of the vacuum p with the ap- 
propriate flow profile is necessary, but also the correct 
implementation of the main contribution at high M. The 
four pion annihilation scenario has to predict an access 
for the high M spectrum that is different from what we 
find find in our temperature dominated evolution. 
Since low temperature (large flow) and high temperature 
(insignificant flow) alter the px spectrum with a differ- 
ent M dependence, both mechanisms cannot lead to the 
same lineshape for different px cuts, hence experimental 
discrimination is possible. 



IV. CONCLUSIONS AND OUTLOOK 

We have discussed a dynamical evolution model for the 
description of 158 AGeV In-In collisions and used it for 
the calculation of the dimuon spectrum measured by 
NA60. 

We have shown that the model is in very good agree- 
ment with the integrated px NA60 spectrum if a 
spectral function including broadening effects and mo- 
mentum dependence based on a ^-functional approach 
is applied. This not only provides a test for the evolu- 
tion model, but also indicates that the vector spectral 
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function found in nature must have similar features, i. 
c. similar px dependence, no strong mass shift and a 
moderate broadening of the p meson. 
We do not claim to have identified the full physics of 
in-medium modifications of the p-meson with this com- 
parison since a scheme that incorporates chiral symme- 
try, self-consistent non-perturbative broadening effects, 
the full three-momentum dependence, finite temperature 
and density as well as the influence of heavier hadrons 
has not yet been achieved. 

We furthermore provided calculations not only of the px 
integrated spectrum, but also of high-p^ and 1ow-pt cuts. 
We expect these data to be published by the NA60 col- 
laboration soon. 

The low-pT analysis allows also a comparison with model 
calculation that do not include a full momentum depen- 
dence of the spectral function. We calculated this spec- 
trum therefore not only for the ^-derivable approach but 
also a thermal nuclear matter calculation [j| and for a 
p-m eson spectral function in a chiral model calculation 
[l8j .The last two models also include further broaden- 
ing via scattering off by baryons such an effect should 
be reflected in an increase of the low M < 770 MeV 
contribution. Substantial broadening is also predicted in 
QCD sum rule approaches with small or no downward 
mass shifts j2j| • The chiral model also includes contribu- 
tions from a broadened in-medium uj and (f>. Considerable 
broadening of the c/>-meson should therefore not only be 
reflected in integrated px, but should also be significant 
at low-pT- 

Discussing the importance of the evolution model we 
found that the px integrated (even more the high-p^) 
spectrum has significant contributions not only from in- 
medium vector p contributions but also from the vacuum 
p. This is especially important in the intermediate mass 
region M « 0.6 — 0.9 GeV where the vacuum rho has 
most of its spectral strength. 

We also found that a purely hadronic evolution is un- 
able to reproduce the data above lGeV invariant mass. 
This was a priori unexpected and is an interesting hint 
that even in the rather small In-In system some amount 
of partonic evolution may contribute to the observed ra- 
diation. Clearly, the contribution is consistent with the 
radiation expected from a QGP, but without further ev- 
idence it would be premature to claim that a QGP has 
been observed in In-In collisions. Data which could help 



to pin down the hadronic freeze-out state (and hence the 
evolution) would clearly help to establish this conclusion 
more firmly. 

It was show that this explanation can be discriminated 
experimentally from a scenario suggested in fl3l | which 
tries to explain the high mass excess by four pion anni- 
hilations. The authors there have shown that they can 
describe the higher mass M > 0.9 MeV part of the in- 
tegrated pt data in their evolution model if they use an 
upper bound for the contribution of four pion annihila- 
tions. It was also shown that considerable mixing has 
to be present even at comparatively low temperatures in 
order to account for the higher mass excess in this sce- 
nario. 

We argue that this explanation can be discriminated from 
our explanation were those higher mass contributions to 
the spectra are mainly due to the early and hot evolution 
stages by considering lincshape in M at different px cuts. 
We provided calculations for the high and low pT-cuts in 
our model. If the apparent good agreement with the 
shape of the preliminary data shown at Quark Matter 
2005 |20j is confirmed, a four pion annihilation scenario 
as the dominant source of the high mass excess can be 
ruled out experimentally. 

We also want to draw the reader's attention to an earlier 
study of us |24| where we argue that an additional exper- 
iment measurement of low mass di-muon production in 
Pb-Pb collisions at 158 AGeV with the same resolution 
as provided by NA60 for the In-In system can put further 
constrains on the theoretical description of the spectra. 
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